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Abstract 

A zigzag boundary between a d^2_y2 and an s-wave superconductor is believed to behave like 
a long Josephson junction with alternating sections of and vr symmetry. We calculate the field- 
dependent critical current of such a junction, using a simple model. The calculation involves 
discretizing the partial differential equation for the phase difference across a long O-vr junction. In 
this form, the equations describe a hybrid ladder of inductively coupled small and vr resistively 
and capacitively shunted Josephson junctions (RCSJ's). The calculated critical critical current 
density Jc{Ha) is maximum at non-zero applied magnetic field Ha, and depends strongly on the 
ratio of Josephson penetration depth Aj to facet length Lf. IfAj/Lj^l and the number of 
facets is large, there is a broad range of Ha where Jc{Ha) is less than 2% of the maximum critical 
current density of a long junction. All of these features are in qualitative agreement with recent 
experiments. In the limit Aj/L/ — > oo, our model reduces to a previously-obtained analytical 
superposition result for Jc{Ha). In the same limit, we also obtain an analytical expression for the 
effective field-dependent quality factor Qj(Ha), finding that Qj{Ha) oc y/jjjia}. We suggest that 
measuring the field-dependence of Qj{Ha) would provide further evidence that this RCSJ model 
applies to a long O-vr junction between a d-wave and an s-wave superconductor. 

PACS numbers: 74.50. +r,74.81.Fa,74.20.Rp 
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I. INTRODUCTION 



Most hole-doped cuprate-based high-temperature superconductors are beheved to be 
characterized by a dj.2_y2 order parameter. Such an order parameter leads to many mea- 
surable consequences, but perhaps the most dramatic are due to so-called vr junctions be- 
tween a dx2-y2 superconductor and either a conventional, s-wave superconductor or another 
dx2-y2 superconductor. These vr junctions can be detected by phase-sensitive symmetry 
measurement spj]. In a vr junction, the Josephson current I and the gauge-invariant phase 
difference 7 across the junction are related by J = /csin(7 -|- 7r)[2!], where /c > is the 
junction critical current. By contrast, a conventional junction satisfies I = IcSin'y. When 
a Josephson loop is formed by a suitable arrangement of d-wave and s-wave superconduc- 
tors, the current-voltage (IV) characteristics differ dramatically from those found when only 
s-wave superconductors are involved. Such phase-sensitive experiments have provided some 
of the most persuasive evidence that the high-Tc superconductors indeed have a dx'2-y2 
superconductor order parameter. 

Recently, the d-wave nature of the order parameter was further confirmed by measure- 
ment of the IV characteristics of long "zigzag" junctions connecting two superconductors: 
a cuprate superconductor with an order parameter of d^^-y^ symmetry, and Nb, an s- 
wave superconductor!^, 4|. In the experiments of Ref. s], the cuprate superconductor was 
YBaaCugOr^^ (YBCO), while in Ref.U Nd2_,;Ce^Cu04-j, (NCCO) was used. The geometry 
of the zigzag junction is shown in Fig. 1 of Ref. y, which also shows the expected orientation 
of the order parameter lobes on each side of the junction. In this configuration, the zigzag 
junction is expected to consist of alternating sections of and vr junctions, as shown in that 
figure. In the experiments, the IV characteristics of this zigzag junction were measured in 
the presence of an external magnetic field applied parallel to the facets of the superconduc- 
tor, e. perpendicular to the plane of the figure. The response of the junctions was, indeed, 
found to be consistent with the expected d^^-y^ symmetry of YBCO and NCCO. 

In this paper, we use a simple numerical model to calculate the critical current density 
Jc{Ha) of a long O-vr Josephson junction as a function of the applied magnetic field Ha- One 
of our goals is to obtain a general picture of how Jc{Ha) behaves in this model as a function 
of the variables \.}/ Lf and Nf, where Lf is the length of one facet of the junction Aj is the 
Josephson penetration depth, and Nf is the number of facets. A second goal is to see how 
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well our simple model agrees with the experimental resultsj^, 4] for O-vr zigzag junctions. A 
final goal is to connect our approach to other treatments of long O-vr junctions. We find that, 
in the limit Aj/Lj ^ 1, the model reduces to a generalized superposition approximation 
previously used to treat such junctions^, Q]. However, even in that limit, we obtain new 
information, namely, an analytical expression for the field-dependent quality factor Qj{Ha) 
of the long junction. 

A characteristic feature of a one-dimensional array of alternating and tt junctions is the 
occurrence of half-integer fiux quanta, in the region where the and vr junctions join. Many 
features of these half-integer fiux quanta have been studied analytically and numerically by 
Goldobin et al. js, 7|, in zero applied magnetic field. Our approach is basically a discretized 
version of this model, but extended to a finite magnetic field. In particular, our model 
does contain the half-quanta which are responsible for some of the observed experimental 
features. 

Qualitatively, our calculated IV characteristics agree fairly well with experiment. For 
example, the measured field-dependent critical current Ic{Ha) of the zigzag junction found 
to be symmetric in Ha, with a principal maximum at nonzero Ha- This behavior is in 
contrast to a long (or long vr) junction, which will have a maximum Ic{Ha) at Ha = 0. Our 
model reproduces these features. We also obtain the result that, when Xj/Lf ^ 1 and the 
long junction has many facets, there is a broad range of Ha where Jc{Ha) is only ~ 1-2% of 
the maximum critical current density. In the experiments of js, J], it is found that Jc{Ha) 
is very small over a wide range of Ha- 

The remainder of this paper is arranged as follows. In Section II, we describe our simple 
model. The numerical results obtained from this model are given in Section III. In Section 
IV, we show that the model reduces to a generalized superposition approximation in the 
limit Xj/Lf — > oo, and we obtain the effective quality factor Qj{Ha) in the same limit. 
Finally, we give a concluding discussion in Section V. 



II. FORMALISM 



We consider a long Josephson junction consisting of alternating and tt sections, or 
facets, as in the zigzag junctions studied experimentally. Let the length of one such facet be 
Lf, and let there be Nf such facets, so that the total junction length is L = NfLj. We use 
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coordinates such that the alternating series of facets runs along the x axis, and the plates 
of the junction are normal to the y axis. In addition we assume that there is a uniform 
magnetic field = i^aZ applied parallel to the junction plates. 

The gauge-invariant phase 7(x, t) across the junction satisfies the following partial differ- 
ential equation: 



the subscripts x and t represent derivatives with respect to position x and a dimensionless 
time ujpt, X being the coordinate along the junction. Aj is the Josephson penetration depth 
and Up is the Josephson plasma frequency, J is the driving current per unit length across the 
junction (assumed uniform), and Jco is the critical current per unit length of a junction at 
zero apphed field {Ha = 0). In terms of the junction parameters, Up = 2eJco/(^C) where C 
is the junction capacitance per unit length. Qj = UpC/G is the junction quality factor, and 
G = is the normal shunt conductance of the junction per unit length. For a junction 
connecting two non-identical materials Xj^ = 87re(Ai + A2 + d)Jc/{hc^)[2l, where Jc is the 
critical current per unit area and Ai and A2 are the penetration depths of the materials on 
either side of the junction. Finally, n{x) indicates whether the point x lies within a or a 
vr facet: n{x) = in a facet and n{x) = 1 in a vr facet. 

In this formulation, the applied magnetic field does not appear explicitly in the differential 
equation; instead, it appears as a boundary condition: 'y{L,t) — 7(0, t) = 27[HadesL/^o + 
[n{L) — n{0)]7[, where $0 = hc/{2e) is the flux quantum. 

The gauge-invariant phase can also be expressed as 



where is the vector potential corresponding to the apphed field H^, and the integral runs 
between the plates of the junction in the y direction. It is convenient to work in a gauge 
such that Aa = HaXy, so that (27r/$o) / A-dl = 27!-Hadesx/^o. d^s is the effective thickness 
of the junction, which we may express in terms of Ai, A2, and the actual Josephson barrier 
thickness d as dcs = d + Xi + X2. 

The equation of motion for the new variable (f){x,t) across the junction may now be 
written as 




Eq. ([1]) is equivalent to eq. (1) of Ref. [tI in the case when Ha is uniform. In this equation. 






(3) 
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with the boundary conditions 0x(O,t) = 4>x{L,t) = 0. For numerical purposes, eq. ([3]) 
is generally more convenient than eq. ([1]), because, in the latter, n{x) is only piecewise 
continuous, and hence n^x is the sum of derivatives of Dirac delta functions. Note also that 
Ha appears explicitly in eq. Q, whereas in eq. ([T]) it appears only as a boundary condition. 

In order to treat eq. numerically, it is convenient to discretize it in space, by breaking 
the long junction into small segments, or "mini- junctions," each of length A {L = NA). 
We also assume that each or tt facet is divided into an integer number Np of such mini- 
junctions, so that Lf = NpA and thus NpNf = N. Denoting the phase difference across the 
i^^ such junction by (pi, we obtain a collection of coupled ordinary differential equations. 
Except for i = 1 and i = N, these may be written 

J 1 

4>i = - sin(0i + 2m f/Np + mTi) + -r^(0i+i - 20^ + + — - — 0i. (4) 

Here we have introduced the frustration / = HaLfdes/^o, defined as the flux per facet in 
units of $0- For i = 1 and i = N, the second term on the right hand side of eq. (jl]) must be 
modified, and the equation of motion becomes 

J 1 • 

0i = - sm{(f)i + 2TTif/Np + TiiTi) + -r|(0i+i - (pi) + —(pi - -p^(t>h (5) 

^ JcO VJ 

and 

)2 J 1 ■ 

(pi = - sm{(pi + 2 TTi f /Np + rii-n) + -A{(pi-i - (pi) + -^(pi - -pz-(pi (6) 

^ JcO hSJ 

for i = 1 and i = N respectively. 

Eqs. (Hj) - ([6]) describe a hybrid Josephson ladder, consisting of N mini-j unctions in- 
ductively coupled together in parallel, and driven by an applied uniform dc current. The 
geometry is shown schematically in Fig. 1. The driving current is in the y direction and the 
field is in the z direction. We use the phrase "hybrid ladder" to distinguish this system from 
the more conventional Josephson ladders, which have junctions on both the rungs and the 



edges js]. The limit of physical interest, however, is the continuum limit, i. e., Np ^ oo. 

III. NUMERICAL RESULTS 

Before describing our results, we make some qualitative remarks about the relation be- 
tween our model and the zigzag junction used in experiments!^, 4|. Although our O-vr 
junction is straight, it does have one crucial feature in common with the zigzag junction, 



namely, intrinsic frustration, even at zero applied magnetic field. We define the frustration 
fi of the plaquette lying between the i*^ and {i + 1)*'^ mini-junction as 

if both junctions are either junctions or vr junctions, and 

= ^ + i (8) 

if one is a junction and the other is a vr junction. If fi is an integer, the energies of all 
the bonds making up the plaquette can be simultaneously minimized, and the plaquette is 
unfrustrated. If fi is non-integer, the bond energies cannot all be simultaneously minimized, 
and the actual bond configuration is a compromise among these; in this case, the plaquette 
is frustrated. 

The effects of frustration are easily understood if A'p = 1. Then the long O-vr junction 
is represented by a hybrid ladder having alternating sections of and vr junctions. Thus, 
according to the above definition, the ladder is frustrated {fi ^ 0) at all fields, including 
Ha = 0, with the exception of / = 1/2, i. e., one-half quantum of flux per facet. At this 
field, the ladder is unfrustrated, and we expect that the critical current per facet will be 
maximum. Moreover, that maximum (per facet) should equal the critical current of a single 
facet at Ha = 0, namely LfJcQ. At all other fields, the critical current will be smaller than 
that maximum value. 

If Np > 1, the ladder has some residual frustration even at / = 1/2. For example, at 
Np = 2 and / = 1/2, each plaquette has fi = 3/4. Thus, although we expect the critical 
current still to be maximum at / = 1/2, that maximum will be smaller than LfJco- Since 
the physical system is really in the continuum limit {Np oo), we generally expect a critical 
current which is smaller than LfJ^o even at the optimum value of / = 1/2. This expectation 
is borne out by our calculations below. 

Although the O-vr ladder configuration of Fig. [1] does have the frustration characteristic 
of a zigzag junction, it still differs geometrically from the one shown in Fig. 1 of Ref. y, 
even in the continuum limit, because the zigzag junction is not straight. The bends in the 
zigzag junction will introduce some pinning which is absent from our model, and therefore 
will probably lead to slightly different IV characteristics than those produced by the model 
of Fig. [H possibly more closely resembling experiment. 



We turn now to our calculations. We have solved eqs. (jl]) - ([6]) for O-vr junctions made 
of various numbers of facets Nf, numbers Np of mini-junctions per facet, facet lengths Lf, 
Josephson penetration depths Aj, and applied fields /. In all cases, we obtain the solutions 
to the coupled differential equations using a standard fourth order Runge-Kutta technique. 
Our goal is to obtain the critical current of a long O-vr junction as a function of these 
parameters, as well as of the Josephson penetration depth Aj, In principle, we are interested 
in the continuum limit, corresponding to A^^p oo. Any effects arising from finite Np will 
be absent from a real O-vr junction, and are therefore artifacts of the finite discretization 
length A. In general, we expect that will vary only on a length scale of Aj. Therefore, 
if Lf <^ A J, we should not need a large Np to approach this continuum limit. Our numerical 
results support this expectation, as we discuss below. 

We first consider a single junction in an external field /. If Aj > L, then the critical 
current per unit length, Jdf), should exhibit a Fraunhofer-like pattern[2]- Fig- 2 of Ref. 
^, which shows their measured Ic{Ha) for a single long junction, clearly exhibits this 
Fraunhofer pattern. 

To calculate this pattern in our model, we apply a specified field / to the junction, then 
ramp up the external driving current density J until a non-zero voltage is recorded, giving 
Jc{f)- f is then increased and theprocess is repeated. Since the value of Qj is not explicitly 
given in the experimental papers 3|, we arbitrarily choose Qj = 6 in all our calculations. 
This value will lead to hysteretic IV characteristics for at least some /, as in the experiments. 

In Fig. [2], we show our calculated Jc{f)/ Jco for a hybrid Josephson ladder of 120 rungs 
with all junctions, open boundary conditions, a penetration depth Aj = 1.3L, and Qj = 6. 
Fig. 0=how= that ou. ladder conBgutation behaves sMariy to the long junction studied 
experimentally (Fig. 2 of Ref. |3|]), and hence, our discretization of the long junction into 
120 mini-junctions is on a sufficiently fine scale to produce the Fraunhofer pattern seen 
experimentally. Comparing Fig. [2] to the experimental Figure, we see that / ~ 1.44 is 
equivalent to Ha ~ l.lyuT. Since all our IV characteristics are obtained on the increasing 
current branch, we have not checked for the expected hysteresis seen in the experiment. 

According to conventional theory[2^, the Fraunhofer pattern in a long Josephson junction 
results from simple superposition: Jc{Ha) = Jco\ sin(7r$/$o)/i,7r$/$o)|; where $ is the fiux 
through the junction area L{d + 2A). This result is obtained 2| by assuming that the phase 
(f){x) [eq. (2)] is independent of x, or equivalently, that (pi is independent of i [eqs. (4) - 
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(6)]. Although, for this particular Xj/Lf, we obtain numerical results very close to this 
asymptotic form, our approach is more general than simple superposition. In particular, 
we do not assume that the 0j's are all the same. Instead, the 0j's, at any applied current 
density J, are always the values which satisfy the coupled differential equations. Thus, in 
principle, our approach gives a different result from superposition, even for this simple case. 
This point is discussed further below. 

In Figs. [31 m and O we show the calculated Jc(/) (obtained on the increasing current 
branch) for a ladder consisting of eight alternating and vr facets, open boundary conditions, 
Qj = 6, and \ j = 2.6L/, but with each facet divided into one, five, and ten mini-junctions, 
respectively. In all three cases, our calculations show Jdf) to be symmetric around / = 0, 
with maxima at / = ±1/2. In the case A = Lf, Jc{f = ±1/2) equals the zero-field critical 
current of a pure junction, but for A = Lf/5 and A = L//10. Jc(±l/2) is smaller than 
this value. This behavior can be understood as follows. For A = Lf, the eight-facet O-vr 
ladder becomes perfectly unfrustrated at / = ±1/2, and thus should behave like the / = 
ladder of junctions. However, in the other two cases, there is still some residual frustration 
at / = ±1/2, and hence Jc(±l/2) is reduced. 

Jc(±l/2) has approximately the same magnitude in Figs. H] and O suggesting that the 
Np = 5 case is already close to the asymptotic limit (A 0), at least for / = ±1/2. In 
general, Jdf) has strong peaks at / = 1/2 ± m, where m is an integer, but these peaks are 
unequal except for A'p = 1. From Figs. H] and [SI we see that the heights of the first three 
peaks in /c(/) do not change when we go from Np = 5 to Np = 10. This is once again a 
consequence of the fact that Np = 10 or even A^p = 5 is close to the continuum limit for 
these values of / and Nf. The model calculations of Figs. [Hand [5] are clearly more realistic 
than that with Np = 1, because they allow </> to vary within a single facet. 

In Fig. [31 but not Figs. [lland[5l Jdf) is periodic with period unity. This periodicity is 
an artifact; it occurs because Nf = 1 in Fig. [HI which causes the equations of motion to be 
periodic with period unity. By contrast, this periodicity is absent if Nf > 1. If Lf = NpA, 
the equations of motion are periodic in / with period Np. But since the physically relevant 
limit is Ap — > oo or A ^ 0, this periodicity is also an artifact of the discretization. 

Another notable feature of Figs. [3E is that each shows exactly six secondary maxima 
between the two principal maxima in Jdf)- These maxima occur at about the same values 
of / no matter how finely the facet is discretized. Moreover, they do not start to overlap 
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as the number of mini-junctions is increased. The maxima can be understood from a gen- 
erahzation of the Fraunhofer pattern described above to the case of Nf alternating and 
TT junctions. This generahzation has been given in eq. (3) of Ref. ^, and, using a shghtly 
different derivation, in Ref. This generahzed Fraunhofer pattern (not shown in the Figure) 
agrees remarkably well with the numerically calculated Jdf) our eight-facet junction. We 
believe the slight deviation is due to the fact that the phases within each facet, in our dy- 
namical model, are obtained directly by solving the equations of motion for the long junction 
rather than by a superposition argument. 

Finally, we comment on another feature of Figs. H] and [5l but not Fig. [3], namely the 
existence of an approximate, but not exact, plateau in Jdf)- In our numerical studies 
for these two ladder systems, we find that the smallest critical current density Jdf) is in 
the range of 1% — 2% of the critical current density Jco of a single facet. Moreover, Jdf) 
remains in this range over a broad range of /. This range of values of 1% — 2% remains 
the same regardless of how small we make the incremental increases in the driving current. 
We believe that this approximate plateau can also be understood from the superposition 
argument mentioned above. More details of this approximate plateau are discussed below. 

To further investigate the behavior of Jdf) in a different range of the ratio Xj/Lj, we have 
carried out additional calculations for a ladder of alternating and vr facets with Qj = 6, 
Nf = 80, and a ratio Xj/Lf = 13, the latter two parameters being the same as in recent 
experiments jsl. We consider three cases: Np = 1, Np = 2, and Np = 4. The results of these 
calculations, shown in Figs. El Wl^a.) and (b), and [HI are quite different from those shown 
previously. Now, for almost all values of / between —0.5 and 0.5, Jdf) ~ 0-01 ^ 0.02Jco- 
Very close to / = ±1/2, Jdf) increases slightly, and exactly at / = ±1/2, Jdf)/Jco = 1 
in Fig. [6l The behavior at / = ±1/2 is consistent with expectations, since, for A'p = 1 the 
ladder is unfrustrated at these values of /. If we consider Np = 2 and Np = 4, but keep 
Xj/Lf = 13, the positions of the peaks remain unchanged and their heights change in a 
predictable manner similar to that already found in Figs. [3H5l In both Figures, as in Fig. 
El there are, as expected, strong peaks in Jdf) at / = 1/2 + m, with m an integer. When 
Np> 2, the peak heights are smaller than unity, as already expected from the results shown 
in Figs, m and El The close resemblance between Np = 2 and Np = 4 suggests that the 
continuum limit is already approached by Np = 4. 

We have carried out many other calculations (not shown here) for other values of Xj/A. 
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In general, we find that this approximate plateau occurs over the widest range of / when 
Aj/A ^ 1 and when Nf ^ 1. Although the experiments are not sensitive enough to 
measure Jdf) in the range of 0.01-0.02Jco, they do show that Jdf) ~ over a broad range 
off. 

We have done various numerical checks to confirm that the plateau is not a numerical 
artifact. These checks are carried out primarily for the results of Fig. 7. To carry out 
the checks, we have typically, for a given /, ramped up J in units of O.OOlJco, where Jco 
is the critical current density of a single facet, rather than O.OlJco as in the calculations 
shown in the other Figures. To compute Jc{f), after ramping up J, we iterated the coupled 
Josephson equations for a dimensionless time Upt = 300, then averaged the voltage over 
the next interval of 200u!pt. Jdf) was taken as the value of J for which the time-averaged 
voltage, as determined in this way, first becomes non-zero. We then incremented / by 0.01, 
reset J to zero, and again incremented the current density in steps of O.OOlJco to find the 
next Jc{f)- In this way, we found that the Jdf) in the plateau region varies between about 
0.012Jco and 0.0167^0- 

The results of these calculations, for a ladder of 80 facets with two mini-junctions per 
facet, are shown in Fig. [7] (a) and (b). Here, we see that Jdf), though small over a wide 
range of /, is not actually perfectly flat, as mentioned above, but instead has small-amplitude 
oscillations as a function of /. Note that, because of the amount of computer time needed 
to do these simulations in which the current is incremented by only O.OOlJco rather than 
O.OlJco, we have done the small-increment calculations only in Figs. 7(a) and 7(b), and not 
Figs. 6 and 8. As in all the other Figures showing Jc{f)/Jco versus /, Jdf) is taken to 
be that value of J above which the time-averaged voltage (V^) jumps up to its value on the 
resistive branch of the hysteretic current- voltage characteristic, i. e. {V) = Jco/G, where G 
is the shunt conductance per unit length. Clearly, Jdf) is indeed non-zero on the plateau, 
and remains nonzero over a broad range of / whenever Nf is large enough. 

We have also calculated Jdf) using the simple superposition theory of Refs. Q and fl. 
The results, shown in Fig. [9], are remarkably similar to those in Fig. [HI Thus, although our 
calculations are based on a numerical solution of the coupled Josephson equations appro- 
priate to this geometry, the field-dependent critical current closely resembles that obtained 
by a simple superposition theory. In particular, the heights of the first two sharp peaks in 
Jdf) calculated numerically, for / > 0, are similar to those obtained using the generalized 
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superposition approximation. 

Finally, to show explicitly how Jc{f)/Jco depends on Aj, we have calculated the critical 
current at zero applied field (/ = 0) for an alternating O-vr ladder with 80 facets, each of 
length one junction [Np = 1), as a function of Aj. The results are shown in Fig. [TUl The 
results of this Figure are calculated incrementing J/ J^o by 0.01 for each value of /. For very 
small Aj the ladder just behaves like many uncoupled small junctions, and Jc(0) = Jco- As 
Aj increases, Jco decreases, until Aj ~ vTO. For Aj > vTO, Jd^) ~ 0.02Jco- Once again, 
this value is not exact, but fiuctuates slightly with Aj. 



IV. GENERALIZED SUPERPOSITION LIMIT 



Since our results closely resemble the simple superposition theory (or "generalized Fraun- 
hofer limit"), we have examined our equations of motion [eqs. (jl]) - (Q] in an effort to un- 
derstand when that model gives the superposition hmit. First, we sum these equations to 
obtain an expression for (1/A^) J2iLi With the definition = (1/A^) J2iLi and the use 
of several trigonometric identities, we finally obtain 

'<i)=-A sin(0 + 6) + J/J,Q - ^/Qj, (9) 

where 

A = ^^C^TS^, (10) 

e = tan-\S/C), (11) 

C = j^Yl cos{2m f/Np + ni7i + 4>i-4>), (12) 

i 

S = sin(27rz//Arp + n^vr + 0^-0). (13) 

i 

The "superposition limit" is appropriate when the phase (pi is independent of i. We 
expect this to be true when Xj ^ Lf, in which case 0, should vary little across the zigzag 
junction. In this limit, 0j — = 0, both C and S are independent of 0, and eq. ([9]) just 
becomes the equation of motion for a single RCSJ with a field-dependent critical current 
and quality factor. Both of these can be obtained with the change of variables ip = (p + 6, 
in terms of which eq. becomes 

= -Asm{^p) + J/J,o - i^/Qj- (14) 
11 



With the further change of variables r' = r/y A, this last equation becomes 

^ = - sin V> + J/(AJ,o) - ^/(Vlgj). (15) 
Eq. (fTSl) describes a single Josephson junction with a critical current density 

Jl{f) = AJ^o (16) 

and quality factor 

Qjif) = ^Qj = Qj , (17) 

where the superscript denotes the superposition approximation. 

To obtain the value of A (in the physically relevant continuum limit A — » 0), we rewrite 
C as C = LimAr^oo(^A)~^ Silif^ cos(27ri//A'p + nj7r] or, converting the sum to an integral, 
C = l/iNfLf) Lfdxcos[2TTfx/Lf + n{x)7c]. The integral can be done analytically. S 
can be obtained in analogous fashion. The resulting expressions for and 5*^, and hence 



A = + and J^if) = AJco, can be shown to reduce to that given in eq. (3) of Ref. |j, 
namely (using our variables as defined earlier) 

{Nf Nf 1 ^ 

J2 sin [n{n - 2nf)] [1 - cos(27r/)] - ^ cos [n{7i - 2nf)] sin(27r/) I (18) 
71=1 n=l J 

and 

{Nf Nf 1 ^ 

J2 cos [n{7T - 2nf)] [1 - cos(27r/)] + sin [n{n - 2nf)] sin(27r/) I (19) 
n=l n=l J 

Thus, our approach does reduce, as it should, to simple superposition in the limit when 
the phases (pi are independent of i. We expect this limit to be appropriate when Aj ^ Lf. 
We also obtain an additional piece of information in this limit, namely, the effective value 
of the junction quality factor Qj{f). Our result for this quality factor is not particularly 
surprising. The reason why the quality factor becomes smaller at certain fields is not that 
the local damping changes with field. Rather, it is because the critical current is reduced 
at certain magnetic fields by cancellation between different parts of the junction, leading to 
an increase in the relative dissipation in the junction at these fields. If \j/Lf is finite, the 
superposition approximation is not exact, because (pi will be dependent on i. In this case, 
we should expect deviations from its predictions for Jdf)- These are, indeed, observed in 
our calculations (and shown explicitly in Fig. 10). 

12 



V. DISCUSSION 



We have presented a simple model for field- dependent critical current density Jc{Ha) of 
a long O-TT junction. An important example of such a system is a zigzag junction connecting 
an s-wave superconductor and one with a d^'^-y-i order parameter. We numerically solve for 
the IV characteristics of this junction by discretizing the equation of motion in space and 
time. The model gives, not only Jc{Ha), but also (in principle) the full IV characteristics 
and its hysteretic behavior when the shunt resistivity is large. 

In the limit Aj — ^ oo, Jc{Ha) for our model reduces to that predicted by a generalized 
superposition modeljs], However, our model is more general, because it applies for finite 
Xj/Lf and hence accounts for the variation of the phase along the junction. Furthermore, 
even in the limit Aj oo, the model gives not only the analytical form for Jc{Ha), as 
previously obtained in Refs. 4] and ^, but also a field-dependent effective quality factor 
Qj{Ha). The field-dependence of Qj{Ha) is due to cancellation effects within the long O-vr 
junction as a function of Ha, rather than any field-dependence of the assumed local damping 
within the junction. 

Our numerical results qualitatively resemble the experiments on d-s zigzag junctions j^, 4|. 
As in the experiment, we obtain a strong maximum in Jc{Ha) at a non-zero magnetic field, 
which is comparable in magnitude to its unfrustrated value, Jco- Furthermore, our results, 
like the experimental ones, depend strongly on the ratio Xj/Lf. Finally, for large Xj/Lf 
and large numbers of facets Nf, we we find that Jc{Ha) is very small over a broad range of 
magnetic field and in the range of 1-2% of Jco- The experiments are not sensitive enough to 
measure such a small Jc{Ha) but they also show broad ranges of Ha where Jc{Ha) is very 
small. The origin of this large region of very small Jc{Ha) is undoubtedly the cancellation 
resulting from superimposing the currents from different parts of the long O-vr junction, as 
implied by the superposition model. The experimental results appear to have some fine 
structure not present in our model. We do not know the origin of this fine structure, but 
speculate that it may be due to the bends in the zigzag junction, which are not included in 
the model. 

Finally, we also comment briefly on the applicability of the RCSJ model to a junction 
connecting an s-wave superconductor to one with d^2_y2 symmetry. One possible concern 
is that, in the d-wave superconductor, since the energy gap vanishes at certain points in 
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k-space, it is relatively easy for a current to excite quasiparticles and, therefore, possibly 
heat the superconductor. Therefore, if a fluxon moves through a long O-tt junction, it might 
easily cause localized heating in the junction. If this heating occurs, it would tend to mask 
the RCSJ behavior of the long O-tt junction. One might also ask how accurately the RCSJ 
model would describe the long junction, even without heating. Clearly, the RCSJ model 
is an idealized approximation of the real junction behavior. But the model seems to be 
a reasonable starting point for possible, more refined approaches. Moreover, the RCSJ 
behavior can potentially be tested, in the limit when the Josephson penetration depth Aj 
is large compared to a facet length. Namely, in this limit, the RCSJ model implies that 
the effective junction quality factor Qj{Ha) is proportional to the square root of the field- 
dependent critical current density Jc{Ha). If experiments show this dependence, this would 
represent evidence that the RCSJ model is indeed applicable to a long O-vr junction. 

The present work suggests several questions which could be further studied in experi- 
ments. For example, it would be of interest to study O-vr zigzag junction when Xj/Lf is 
smaller than unity. In this range, we would expect substantial departures from the gen- 
eralized superposition model, as suggested by Fig. [10] for / = 0. In addition, it would be 
desirable to test the predicted dependence of Qj{f) on / experimentally in the superposition 
limit. Since Jc(/) may be small over a broad range of /, eq. (ITTI) suggests that a long O-vr 
junction would behave as if overdamped at these values of / even if it is underdamped and 
hysteretic when Jc{f) is large. It would be of interest to test this hypothesis experimentally. 
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FIG. 1: Schematic view of Josephson ladder used in present calculations. The ladder has open 
boundary conditions; u and £ represent the upper and lower edges of the ladder. A uniform current 
density J per unit length is applied in the y direction. This corresponds to a current / = JA applied 
across each rung in the ladder. An external magnetic field is applied in the z direction. 
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FIG. 2: Calculated critical current per unit length Jc{f) for a junction of length Lf, plotted 
as function of the flux / through the facet. / is given in units of a flux quantum $o = hc/{2e). 
The facet is represented as a hybrid ladder of 120 small junctions, as described in the text. The 
Josephson penetration depth Aj = 1.3L/ and the quality factor Qj = 6.0. Jco is the critical current 
per unit length of a junction at zero applied field. 
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FIG. 3: Critical current per unit length Jdf), in units of Jco, plotted as a function of the magnetic 
flux per facet /, for an eight-facet junction of alternating and tt facets, open boundary conditions, 
A J = 2.6L/ and Qj = 6.0. The junction is modeled as a hybrid Josephson ladder with one mini- 
junction per facet, as described in the text. 



18 




-1 -0.5 0.5 1 



f 

FIG. 4: Same as in Fig. [3l but with five mini-junctions per facet. 
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FIG. 5: Same as in Fig. [3l but with ten mini-junctions per facet. 
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FIG. 6: Same as in Fig. [3l but with 80 facets, one mini-junction per facet, and Aj = 13L 
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FIG. 7: (a) (top): Same as in Fig. [U but with two mini-junctions per facet. In contrast to Fig. [U 
we have incremented the current density in increments of O.OOlJco to calculate the critical current. 
With this fine scale, small oscillations in Jc(/) are clearly apparent in the "plateau" region, which 
is not perfectly flat, (b) (bottom): same as (a), except that we have magnified the vertical scale 
to make the oscillations more visible. 22 
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FIG. 8: Same as in Fig. [U but with four mini-junctions per facet. As in Fig. [6l we increment the 
current density in intervals of O.OlJco- At this resolution, most of the plateau region looks flat, but 
we expect that, if calculated with increments of O.OOlJco it would show oscillations as in Fig, [71 
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FIG. 9: Jc{f)/JcO as obtained from the superposition approximation of Refs. [J] and [5| for a lon^ 
O-vr junction having 80 facets 
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FIG. 10: Critical current per unit length Jc{f), in units of Jco, for / = and an 80-facet O-tt 
junction modeled as in the text for / = 0, assuming one mini-junction per facet, and plotted as a 
function of Xj/Lf. 
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